clear;clc;
close all;
FimageQ=load('Fimage_Qacoustic.mat');

FimageA=load( 'Fimage_acoustic.mat');

FimageA_Qdata=load( 'Fimage_acoustic_attenuated.mat');

[zz,xx]=size(FimageQ.Fimage);
FimageQ.image=zeros(size(FimageQ.Fimage,1),size(FimageQ.Fimage,2));
FimageQ.image(2:end,:)=FimageQ.Fimage(2:end,:)-FimageQ.Fimage(1:end-1,:);

dx=8;
dz=8;
x_tick=[0:100:xx];x_label=dx*x_tick;
z_tick=[0:100:zz];z_label=dz*z_tick;
Font=14;

plotimage(FimageA.Fimage);
title('acoustic');
set(gcf,'units','normalized','position',[0.1 0.1 0.6 0.6]);
axis([0,xx,0,zz]);set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('x (m)','FontSize',Font);
ylabel('z (m)','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'XTick', x_tick);
set(gca,'XTickLabel',x_label,'FontSize',Font);
set(gca, 'YTick', z_tick,'FontSize',Font);
set(gca,'YTickLabel',z_label);

plotimage(FimageA_Qdata.Fimage);title('acoustic with attenuated data');
set(gcf,'units','normalized','position',[0.1 0.1 0.6 0.6]);
axis([0,xx,0,zz]);set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('x (m)','FontSize',Font);
ylabel('z (m)','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'XTick', x_tick);
set(gca,'XTickLabel',x_label,'FontSize',Font);
set(gca, 'YTick', z_tick,'FontSize',Font);
set(gca,'YTickLabel',z_label);

plotimage(-FimageQ.image);title('Q acoustic');
set(gcf,'units','normalized','position',[0.1 0.1 0.6 0.6]);
axis([0,xx,0,zz]);set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('x (m)','FontSize',Font);
ylabel('z (m)','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'XTick', x_tick);
set(gca,'XTickLabel',x_label,'FontSize',Font);
set(gca, 'YTick', z_tick,'FontSize',Font);
set(gca,'YTickLabel',z_label);



%%%%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%%%
%%%            plot trace             %%%%%%%%%%
%%%%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%%%
xi=300;
figure;set(gcf,'units','normalized','position',[0.1 0.1 0.6 0.4]);
plot(FimageA.Fimage(:,xi),'r','LineWidth',1);
hold on;
plot(-FimageQ.image(:,xi)*0.5,'--b','LineWidth',1);
plot(FimageA_Qdata.Fimage(:,xi),'k','LineWidth',1);
xlim([0,zz]);%set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('z (m)','FontSize',Font);
ylabel('Amplitude','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'xTick', z_tick,'FontSize',Font);
set(gca,'xTickLabel',z_label);

xi=600;
figure;set(gcf,'units','normalized','position',[0.1 0.1 0.6 0.4]);
plot(FimageA.Fimage(:,xi),'r','LineWidth',1);
hold on;
plot(-FimageQ.image(:,xi)*0.75,'--b','LineWidth',1);
plot(FimageA_Qdata.Fimage(:,xi),'k','LineWidth',1);
xlim([0,zz]);%set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('z (m)','FontSize',Font);
ylabel('Amplitude','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'xTick', z_tick,'FontSize',Font);
set(gca,'xTickLabel',z_label);

xi=800;
figure;set(gcf,'units','normalized','position',[0.1 0.1 0.6 0.4]);
plot(FimageA.Fimage(:,xi),'r','LineWidth',1);
hold on;
plot(-FimageQ.image(:,xi)*0.75,'--b','LineWidth',1);
plot(FimageA_Qdata.Fimage(:,xi),'k','LineWidth',1);
xlim([0,zz]);%set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('z (m)','FontSize',Font);
ylabel('Amplitude','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'xTick', z_tick,'FontSize',Font);
set(gca,'xTickLabel',z_label);
